Long-distance entanglement and quantum teleportation in XX spin chains 
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Isotropic XX models of one-dimensional spin- 1/2 chains are investigated with the aim to elucidate 
the formal structure and the physical properties that allow these systems to act as channels for 
long-distance, high-fidelity quantum teleportation. We introduce two types of models: I) open, 
dimerized XX chains, and II) open XX chains with small end bonds. For both models we obtain 
the exact expressions for the end-to-end correlations and the scaling of the energy gap with the 
length of the chain. We determine the end-to-end concurrence and show that model I) supports 
true long-distance entanglement at zero temperature, while model II) supports "quasi long-distance" 
entanglement that slowly falls off with the size of the chain. Due to the different scalings of the 
gaps, respectively exponential for model I) and algebraic in model II), we demonstrate that the 
latter allows for efficient qubit teleportation with high fidelity in sufficiently long chains even at 
moderately low temperatures. 

PACS numbers: 03.67.Hk, 03.67.Mn, 75.10.Pq 



I. INTRODUCTION 



The crucial role of entanglement as a fundamental re- 
source for quantum information tasks that transgress the 
classical limits is particularly evident and has been ex- 
perimentally verified in such protocols like teleportation 
IE j cryptography and secure key distribution [l| , and 
quantum communication Typically, the "most natu- 
ral way" to create entanglement between parties is by 
means of direct interactions. Since, intuitively, large 
amounts of entanglement should be associated to the 
presence of strong correlations, low-dimensional systems, 
as, for instance, spin chains, offer a natural source of en- 
tanglement. However, in most systems with short-range 
interactions, the entanglement between a pair of particles 
decays rapidly with the distance. For instance, in the 
Ising model with transverse field || the two-spin concur- 
rence vanishes for distances larger than two neighboring 
sites, while in the Heisenberg model @, it is restricted 
only to single nearest neighbors. A first exception to this 
behavior was found by Amico et al Q, who established 
that close to factorization points, the range of entangle- 
ment grows indefinitely. However, in this case, the en- 
tanglement strength between two spins rapidly vanishes 
as the range increases. In fact, this behavior appears 
to be natural, recalling that monogamy of entanglement 
Q implies that if entanglement is enhanced between two 
given subsystems, then the entanglement between any of 
the two subsystems and a third one must necessarily be 
suppressed. 

From a general quantum informatic perspective, a 
much desired goal would be the ability to create a large 
amount of entanglement between distant subsystems, si- 
multaneously avoiding direct interactions between them 



and single-subsystem addressing. Along this line of 
thought, the first requirement can be satisfied by intro- 
ducing the concept of localizable entanglement, with the 
aim of exploiting spin chains as quantum channels 0]. 
The localizable entanglement measures the average en- 
tanglement that can be concentrated on a pair of dis- 
tant subsystems by performing optimal local measure- 
ments onto the rest of the system. More recently, it was 
shown that the ground state of some spin models with fi- 
nite correlation length defined on one-dimensional chains 
with open ends can support large values of long-distance 
entanglement between the end points of the chain [l(|. 
This approach, i. e. to look for systems whose ground 
state can support long-distance entanglement prima fa- 
cie, without the need of performing operations and mea- 
surements, is clearly very appealing. However, it ap- 
pears that models which exhibit true long-distance en- 
tanglement are characterized by energy gaps above the 
ground state that vanish exponentially as the length of 
the chain is increased [l(|. Hence, this interesting phe- 
nomenon seems doomed to survive only in the physi- 
cally unattainable situation of zero temperature. In the 
effort to overcome this problem, it was shown that a 
kind of long-distance entanglement, very slowly decreas- 
ing with the length of the chain, can be supported by 
the ground state of spin models with infinite correlation 
length defined on one-dimensional open chains with small 
end bond interactions, and that these systems allow for 
robust finite-temperature teleportation across finite dis- 
tances [ll| . This "quasi long-distance entanglement" be- 
tween the end points of the chain is associated to an en- 
ergy gap that vanishes algebraically with the size of the 
system. As a consequence, it is more resilient to thermal 
excitations, and can be used to engineer realistic proto- 
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cols of qubit teleportation, as has been explicitly demon- 
strated, by numerical DMRG simulations, in the case of 
the Heisenberg (XXX) chain with small end interactions 

In the present paper we will show that both true long- 
distance and quasi long-distance entanglement can be 
supported by very simple isotropic XX models of open 
spin chains. These models present some interesting ad- 
vantages over the ones that have been previously studied. 
At variance with the Heisenberg case, they are amenable 
to exact analytical treatment both in the case of alter- 
nating couplings, corresponding to dimcrization of the 
ground state and true long-distance entanglement, and 
in the case of small end bonds, corresponding to a ground 
state that supports quasi long-distance entanglement. In 
this way, one achieves a full grasp of the physical mech- 
anism responsible for long and quasi-long distance en- 
tanglement, and the possibility to identify unambigu- 
ously the optimal range of parameters for teleportation 
with maximal fidelity. Moreover, having in mind pos- 
sible experimental realizations, for instance in suitably 
engineered optical lattices XX chains are in princi- 
ple more easily realizable or simulatable than Heisenberg 
interactions. The exact solvability of the open-end XX 
chains is possible thanks the methods introduced by Lieb 
and coworkers [HI . We will first study the fully dimerized 
XX open spin chain with alternating couplings, discuss 
the exact behavior of the end-to-end concurrence in the 
ground state at zero temperature, and show analytically 
that for this kind of model the energy gap between the 
ground and the lowest excited state falls off exponentially 
with the size of the system. Therefore, this system can- 
not be exploited for quantum teleportation in realistic 
situations at finite temperature. We then study the XX 
open spin chain with small end bonds and determine an- 
alytically the exact expression for the zero temperature 
end-to-end concurrence in some physically relevant lim- 
its. We determine the exact scaling of the energy gap 
with the length of the chain and show that it dies off 
algebraically with the size of the system. Finally, we 
determine the behavior of the teleportation fidelity as a 
function of the temperature for different strengths of the 
small end bonds, and show that qubit teleportation with 
fidelities well above the classical threshold, and in some 
cases close to unity, is supported even at moderately low 
temperatures. 

The paper is organized as follows: In Section II we 
introduce the general XX spin chain Hamiltonian with 
arbitrary site-dependent couplings, and solve it analyt- 
ically for the end-to-end two-point reduced density ma- 
trix, correlation, and concurrence (entanglement of for- 
mation). We calculate the fully entangled fraction and 
establish the analytical expression for the fidelity of tele- 
portation of an unknown qubit state between the two 
end points of the chain. In Section III we specialize the 
model to the case of perfectly alternating couplings and 
fully dimerized ground state. We determine the single 
particle dispersion law, establish the exponential scaling 



behavior of the energy gap, and determine the analytical 
expressions for the ground state end-to-end concurrence 
and fidelity. We show that, depending on the values of 
the coupling strength, high or even maximal entangle- 
ment and fidelity can be achieved for chains of arbitrary 
length. In Section IV we turn to the case of open XX 
chains with small end bonds, establishing the single parti- 
cle dispersion law, the algebraic scaling of the energy gap, 
and the analytic expressions for the concurrence and the 
fidelity. We then investigate the behavior of the end-to- 
end teleportation fidelity as a function of the temperature 
for different values of the small end bonds, and conclude 
that teleportation with unit fidelity at very low temper- 
ature and with high fidelity above the classical threshold 
at moderately low temperature are both supported by 
the XX channel with small end bonds. Finally, in the 
conclusions we summarize our findings and discuss some 
outlooks on possible future developments along this line 
of research. 



II. THE GENERAL MODEL 

As already anticipated in the introduction, we will fo- 
cus our analysis on one-dimensional lattices with open 
ends, described by XX models with different types of 
nearest neighbor interactions. Such models are all spe- 
cial instances of the general XX Hamiltonian 

L-l 

H = Y,Ji(s?sf +1 +s»sy +1 ) , (i) 

i=l 

where Ji is the interaction strength between nearest 
neighboring sites i and i + 1, Sf are the spin opera- 
tors defined at site i, and L is the total number of sites 
(spins) or the length of the chain. The spectrum of this 
Hamiltonian can be determined exactly by a straightfor- 
ward application of the standard methods introduced in 
Ref. [131 ] . The first step in the procedure is to perform a 
Jordan- Wigner transformation [l4| , 

S! = ch (2) 

where Sj = 5J ± iS? . As a result, the Hamiltonian (TTJ 
is mapped in the free fermion Hamiltonian 

h = - ^ (4 c i+i + 4+i c i) = ctMc » ( 3 ) 

i=l 

where = (c\ , . . . , c\\ (c) is the vector of the L creation 
(annihilation) operators, one for each site of the lattice, 
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and the adjacency matrix M reads 

J i ••• 

Ji J 2 



M=- 



J 2 







• ■ J L -2 
JL-2 J L _! 

Ji_i 



(4) 



We will investigate two particular realizations of Hamil- 
tonian {!]) that, as we will show, in the ground state allow 
for high or even maximal entanglement between the end 
spins of the chain, and thus naturally provide a channel 
with high or even unit fidelity for qubit teleportation. 
To evaluate the teleportation fidelity, one needs to deter- 
mine the spin-spin concurrence (entanglement of forma- 
tion) between the end points of the chain. This quantity 
can be computed exactly for any two-qubit state (pure 
or mixed) , thanks to the celebrated formula of Wootters 
[HI], and the task is left to obtain its explicit expression 
in the reduced state of the two end-point spins. To this 
purpose, we need to calculate explicitly all the possible 
forms of two-point correlations in the ground state. The 
Hamiltonian ([T]) is symmetric under rotations of the spins 
around the z-axis, so that the only nonvanishing corre- 
lations are (SfSj) = (SfSlj), (S?S?) and (Sf). In the 
absence of external magnetic fields, 7r-rotations around 
the x and y axes are symmetries of the model, which 
additionally implies (Sf) — at every site. Thanks to 
the aforementioned symmetries, the two-point reduced 
density matrix pij, obtained by tracing the full density 
matrix of the system over all sites except the pair {i,j}, 
has the form 



I 



+ (SfSf) (a x ® a x + a y ® a v ) + (SfS*)a z 



(5) 

where a x,v,z are the Pauli matrices and (•) is the ground 
state average at temperature T = 0, or the thermal one 
at finite temperature f3 — (fc^T) -1 with respect to the 
Gibbs state p = e^Z' 1 . We are interested in the case 
in which i and j are the two end points of the chain. In 
this instance, we have 



SiS L + S 1 s£ — - 
SiS' L = Ucx ~ 



C L CL-~ 



(6) 



where A/" = Yli=i c \ c i IS the total number operator. Us- 
ing Wick's theorem and taking into account that (c\ci) — 
1/2, we obtain 

{St SI + S^S+) = -e^/ 2 ((c\c L ) + (4 Cl >) , 
(S!Sl)=-(c\c L )(cl Cl ). (7) 

Setting x = (c\cl) we see that the end-to-end reduced 
density matrix depends uniquely on this parameter. 



All the physical informations about model ([3]) can now 
be obtained by diagonalizing the one-body matrix M . 
Let £fc be the eigenvector with eigenvalue Afc, where k is 
a quasi-momentum label. Then, passing to new fermionic 

operators via the transformation Cj = X^/c^^ c fc> the 
Hamiltonian takes the form 



H = A fc 4cfc 



k 



(8) 



The evaluation of the two-point correlation x is then 
straightforward, and one obtains 



k . q 



(1)AL) 1 



for T > 



(9) 



For reduced states of the form |5]), the end to end con- 
currence is easily computed, and we have fl6j 



C'i _i = 2 max { 0, x 2 + \x\ - - 



(10) 



The above expression of the concurrence is nonvanishing 
for \x\ > (v / 2-l)/2 ~ 0.207, and it reaches the maximum 
value C hL = 1 for \x\ 1/2. 

It is natural to expect that the existence of a strong 
quantum correlation between the two end spins of the 
chain can be conveniently exploited for performing tasks 
in quantum information, in particular considering tele- 
portation schemes. In the standard quantum teleporta- 
tion protocol, two parties A and B share a maximally 
entangled state (Bell state). Party A holds also a third 
qubit, whose unknown state is to be teleported. If the 
two end points of our XX chain share a highly entangled 
state, that in some limit may even be asymptotically close 
to a Bell state, they can be identified as the two parties, 
sender and receiver, for a long-distance, high-fidelity tele- 
portation protocol. The efficiency of a quantum channel 
in teleporting an unknown state is quantified by the fi- 
delity / between the output and the input states, aver- 
aged over all input realizations. The fidelity depends on 
the actual properties of the entangled resource p\^ (Cfr. 
Eq. ([5])) shared by the end spins of the chain. In fact, it 
has been demonstrated that the optimal fidelity depends 
only on the "fully entangled fraction" Ff u u , according to 
the formula / = (2F fua + l)/3 0. The fully entangled 
fraction is defined as the fidelity between the resource 
p± t L and a maximally entangled state, maximized over 
all possible maximally entangled states. For states of 
the form Eq. ([5]) it can be easily computed, and reads 
Ff u ii = \ + \x\ + x 2 [l8[. The associated teleportation 
fidelity is thus 



/ 



2(l + \x\+x 2 )+l 



(11) 
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This expression highlights the crucial interplay between 
entanglement and efficiency in quantum information pro- 
tocols. In fact, due to the high symmetry of states of 
the form Eq. ([5]), a nonvanishing entanglement implies a 
nonclassical teleportation fidelity exceeding the classical 
threshold 2/3, and viceversa. In the limit |x| — ► 1/2 of 
maximally entangled resource, the maximum teleporta- 
tion fidelity reaches unity. 



III. LONG-DISTANCE ENTANGLEMENT 



In analogy with previous work on the dimerized 
Heisenberg model [lOj], we consider first the open end 
XX chain with bonds of alternating strengths (1 — S) 
(weak bond), and (1 + 6) (strong bond), with < 5 < 1. 



L-l 



(12) 



Choosing L even and < 6 < 1 assures that the spins 
at the end of the chain interact with a weak bond of 
strength (1 — 8) with their respective neighbors. 

Let us first comment on the general features of the 
model (fT2"|) in the thermodynamic limit. For (5 = the 
model in the fermionic picture reduces to a simple tight 
binding with dispersion Afc = J cos (fc). The ground state 
is given by a half filled band and vanishingly small ex- 
citations are present near the Fermi points. Upon the 
introduction of a nonvanishing S translational invariancc 
by one site is broken, the Brillouin zone is correspond- 
ingly halved, and a gap of width 6 opens up in the single 
particle spectrum. Two bands develop with dispersion 



A fe>± = ±Jv/cos 2 (k) + 4<5 2 sin 2 (k) 



(13) 



Upon increasing the value of 5, the growing difference 
between the alternating coupling constants allows the 
creation of dimers between pairs of strongly interacting 
spins. The dimers, in turn, are very weakly interacting 
with each other. Hence, because the two end spins of 
the chain interact very weakly with their nearest neigh- 
bors, dimerization and monogamy of entanglement force 
the creation of a strong quantum correlation between the 
end points, in analogy with what happens in the Heisen- 
berg case (lOj . 

The above qualitative picture can be supplemented by 
a detailed quantitative analysis both in the thermody- 
namic limit and at finite length L. The eigenvectors of 
the adjacency matrix M are given by [l|| 



f U) 



A k sin (kj) 
H k A k sin (k (L - j 



j even , 
1)) j odd , 



(14) 



where /ifc = ±1 is the parity (left-right symmetry) of the 
eigenstate and Aj~ is a normalization constant that 
reads 



A k = 2 



2 L + l- 



sin2fc(L + i; 
sin (2k) 



-1/2 



(15) 



The quasimomenta are given by the solution of the fol- 
lowing equation 



sin k (L + 2) 
sin (kL) 



1 



1 

a 



(16) 



For each quasimomentum wc have two eigenvectors cor- 
responding to opposite parity. Hence we seek for L/2 
solutions of HID. When a < L/ (L + 2) (roughly 5 > 0, 
when L is large) we have L/2 — 1 real solutions for which 
the energy dispersion has the form (fT3")l . The missing so- 
lution is complex (and plays a very important role) and 
has the form ko = n /2 + ip. For this mode the equation 
reads 



cosh (2p) + coth (Lp) sinh (2p) - 
At leading order the solution is given by 



o2p 



1 



- (l-a 2 )^- 1 



(17) 



(18) 



The appearance of a complex quasimomentum is crucial 
for the existence of long-distance entanglement (LDE) in 
the ground state and can be understood as follows. The 
eigenvectors £,k ,± are localized at the boundaries of the 
chain. For instance, for site j close to the left border 

of the chain (j ~ 1), one has ^ ± ~ e ( L -J+ 1 )/C, with 

a localization length £ that, for an asymptotically large 
chain, reads 



C = 2/ In [(1 + 6)/ (1-6)] 



(19) 



The energy of this mode, for both value of the parity, is 
exponentially small: 



A fe0i± ~±2(l-a)e 



-L/C 



(20) 



and we thus see that the onset of LDE in the dimerized 
ground state of the XX chain with open ends is strictly 
linked to the existence of an energy gap above the ground 
state that vanishes exponentially fast with the length of 
the chain. 

At zero temperature and zero external field, all nega- 
tive energy modes are filled, up to A ko _. The mode k is 
responsible for the appearance of a state localized at the 
end of the chain in the many-body ground state. This 
means that the ground state resemble ^ ~ ® V'rest- 
As we have dimer order at sites 1 and L, we expect 
l to be highly entangled. Because factorization will 
not be exact, i/ji.l will not be exactly a pure state but 
rather a weakly mixed reduction p\ t L that becomes max- 
imally entangled only in the asymptotic regime of perfect 
dimerization. In fact, knowing the eigenvectors of the ad- 
jacency matrix, we can evaluate exactly the end-to-end 
LDE between the borders of the chain. Using Eq. (Q2 
the fermionic correlation function x reads 



(c\cl) 



= ]T ^At[sin(Lk)Y 



(21) 



k, A fc <0 
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To evaluate this sum we first isolate the contribution 
coming from the complex momentum. The remaining 
set of terms defines a sum with alternating signs of a 
periodic function analytic on the real axis. The rate 
of convergence to its asymptotic vanishing value is dic- 
tated by the width of the largest strip, around the imag- 
inary axis, where the function is analytic. The inverse 
of this width is precisely given by the localization length 
Eq. (|19p . Such contribution is exponentially suppressed, 
and the final result is 



(-l) i/2 i(l-a 2 )+0(e-VC) 



(-1) 



L/2 



26 



O e- L /< 



(22) 



We have checked this approximation against ex- 
act numerical data and the results are plotted in 
Fig. [1] It turns out that the remainder has the form 
(— l) L ^ 2 A a L 2 cxp (— L/C) +•••, where Aq is a positive 
constant. This implies that, in modulus, the asymptotic 
value of the end-to-end correlation is in fact approached 
from below: In this type of systems, correlations increase 
as the length L of the chain is increased. 

We see that, for any nonvanishing value of the dimer- 
ization, the ground state of the system develops a con- 
stant, non-zero correlation between the end points of 
arbitrarily large chains, a phenomenon known as sur- 
face order. Using Eq. (flU)) , we deduce that this sur- 
face order, present whenever 6 > (and 8 < 1) allows 
for LDE between the end spins as soon as 5 > 8q — 

(l - / (l + V% - ~V$\ = 0.132, i.e. for mod- 

erately low values of the dimerization. The correspond- 
ing threshold localization length is Co = 7.479. The lo- 
calization length C rapidly decreases from the threshold 
value Co with increasing 5 > So- As a consequence, the 
asymptotic regime for the corresponding end-to-end con- 
currence (LDE) is reached already for chains with a size 
of few tens of sites. Finally, by substituting in Eq. (jTUJ) 
and in Eq. ifTTj) the expression of x Eq. (|2"2")1 , we obtain 
the analytic expressions of the LDE (end-to-end concur- 
rence) and of the teleportation fidelity as functions of the 
dimerization ratio a = (1 — 5)/ (I + 8) in the asymptotic 
regime of chains of large size L 3> Co- One has: 



2 max < 0, 



1 



a 

T 



/ 



1 



a 

T 



(23) 



(24) 



In Fig. [2] we report the behavior of the LDE between 
the end points of the chain as a function of the dimeriza- 
tion parameter 8 for large chains. We see that as soon 
as the dimerization is above the threshold value 8q, the 
LDE grows very rapidly with 8, reaching saturation in 
the limiting situation <5 — > 1. 

Before ending this section we should remark that, al- 
though extremely interesting, the true LDE picture, de- 
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Figure 1: (Color online) Absolute value of the zero- 
temperature correlation xl = \{c\cl)\ as a function of the 
size L of the chain according to the scaling xl = x x — 
A L 2 exp(-L/C). The different curves reproduce ln[(xoo — 
xl)/L 2 ] as a function of L for different values of the dimeriza- 
tion parameter 8. Symbols denote the exact numerical values. 
The values of £ are obtained from Eq. (|19|) . 




Figure 2: End to end concurrence Ci,l, for long chains of 
length L S> Co in the dimer model (I12|l . as a function of the 
dimerization parameter 8. Nonvanishing LDE is achieved and 
is rapidly growing as soon as 8 > So = 0.132. 



rived at zero temperature, does not survive, even qual- 
itatively, at finite temperature. We have seen that the 
presence of LDE in the dimerized XX is connected to 
the presence of an eigenstate of the adjacency matrix 
with complex quasimomentum. From Eq. (|20|) it follows 
that the energy gap between the ground state and the 
lowest excited state vanishes exponentially as the size of 
the chain increases. Hence, even for small chains and 
very low temperatures, the first excited levels get signifi- 
cantly populated, contributing to the correlation x with 
values of opposite sign with respect to that of the ground 
state. The total effect is to lower x down to zero even 
at very low temperatures. In conclusion, we have deter- 



mined that LDE can be supported at zero temperature by 
relatively simple, dimerized XX chains with open ends. 
However, in close analogy with more complex systems 
analyzed previously [HI ] , the zero temperature LDE does 
not survive as soon as temperature is switched on. 



IV. QUASI LONG-DISTANCE 
ENTANGLEMENT 



In this section we discuss a second type of model that 
can support highly entangled reduced states at the end 
points of a spin chain. In all generality, spin models that 
can allow for strong end-to-end correlations are charac- 
terized by interactions between the end points and their 
nearest neighbors that are smaller compared to the inter- 
actions in the bulk of the chain. Otherwise, if the system 
does not meet this criterion, the end points would be- 
come strongly entangled with their neighbors, excluding, 
due to monogamy constraints Q, the possibility of LDE. 
We then consider a model of open XX spin chain formed 
by L — 2 spins with uniform coupling strengths, plus two 
weakly interacting probes placed at the two end points. 
Such a models is described by the following Hamiltonian 

L-2 

H = jY,[S*S* +1 + SfS! +1 (25) 

i—2 

, \ (qxqx i ay ov i qx qx , qy qy\~\ 



where < A < 1. The presence of the A-term in Eq. 
can be understood as a kind of generalized boundary con- 
dition, reducing to the standard XX model with uniform 
couplings and open ends for A = 0, 1. It is then not 
surprising that the eigenvalues of the corresponding ad- 
jacency matrix M have the form Afe = J cos (A;). The 
eigenvalue equation determining the quasimomenta for a 
generic A reads [2(| 



[i k cot (k) 



L - 1 
cot I — - — k 



A 2 



2- A 2 ' 



(26) 



where ji^ = ±1 is again the parity of the corresponding 
eigenstate. The eigenstate associated to the quasimo- 
mentum k has the following components: 



(1) 



« 



Ak 
1 

+(1 - A 2 ) sin((i - l)fc)) , 1< i < L , 
A 



• sin(fc) , 
(sin((i + l)fe) 



fi— sin(fc) , 



where the normalization Ah reads 



Al = (L-l) 



2(l-A 2 )cos 2 (fc) + ^ 



(27) 



+ 2A 2 -A 4 . (28) 




Figure 3: (Color online) End-to-end concurrence Ci,l at zero 
temperature as a function of A for different lengths L of the 
chain. Curves from right to left are, respectively, for L — 26; 
L = 50; L = 100; L = 200; and L = 400. 



Eq. (|26|) admits L distinct real solutions in the interval 
(0, 7r) for A ^ 0. Using Eq. (|27p the zero temperature 
end-to-end correlation x reads 



E 



7r/2<fe<7T 



A . a 



(29) 



The sum is limited to the L/2 quasimomenta that lie 
in the half-interval ir/2 < k < ir. This condition takes 
into account the fact that, at exactly zero temperature, 
all levels below the Fermi energy are occupied, while all 
the others are empty. The value of the end-to-end con- 
currence Ci t L is readily obtained by inserting expression 
(|29| for the end-to-end correlation x in Eq. ([TO]) . In 
Fig. [3] we show the behavior of the end-to-end entan- 
glement as a function of A, for different lengths of the 
chain. At fixed length L of the chain, C\ l — » 1 as 
A — > 0. At a fixed value of A, the end-to-end concur- 
rence slowly dies off as L — > oo. The XX spin chain 
with small end bonds is thus characterized by what we 
may name "Quasi Long-Distance Entanglement (QLDE) 
because, at variance with the dimerized case considered 
in the previous section, the end-to-end entanglement, al- 
though able to reach asymptotically maximal values in 
the limit of vanishing end-couplings, slowly decreases as 
the size of the chain is increased. These considerations 
can be made quantitative by investigating in detail the 
exact expressions for x that are obtained from Eq. (T2"9"|) 
in the two limiting cases A — > 0, 1. 

In the limit of vanishing small end bonds (A — > 0), the 
quasimomenta assume the expressions 



tan 



k n 



L - 1 



- A 



L-l) 



L - 1 



l,2,...,L-2, 



± 



(30) 



Isolating the dominating contributions in the sum 
and keeping terms up to second order in A, in the limiting 
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interpolating 
approximation X ~ 1 
approximation X ~ 




Figure 4: (Color online) Comparison between the exact nu- 
merical evaluation and the analytic approximations for the 
end-to-end correlation function x = {c^cl) for a chain of 
L — 100 sites. The intepolating curve is given by Xi nt = 
M) V2 / ( 2 + 4cLA 2 ). It is exact up to O (A 2 ) when A -> 0, 
but retains its validity also for higher values of the coupling. 



case A — > 0, one obtains 



x\-^q = (-1) 



L/2 



A i 



2C 

^2~ 



-O (A 4 



, (31) 



where C is Catalan's constant C = 0.915, and the con- 
stant in brackets is c = 0.435. The above result pro- 
vides the quantitative content of the qualitative picture 
sketched above: At fixed size L of the chain, one can 
always choose A small enough so that |x| — * 1/2 and 
the end-to-end concurrence approaches unity. Depending 
on the length of the chain, the condition to be satisfied 
to achieve large values of the end-to-end concurrence is 
A < 1/VI. 

For completeness, let us now consider the opposite sit- 
uation of strong end bonds A — > 1. In this case, all the L 
quasimomenta can be approximated by k n ~ nn/ (L + 1) 
with n — 1, . . . , L. Exploiting the fact that, at leading 
order, the alternating Riemann sum of a function is equal 
to (one-half of) the value of the function at the extrema, 
we obtain 



(-1) 



L/2 



1 



A 2 (L-3) + 4 



(32) 



In Fig. 0] we compare the exact expression of X clS cL func- 
tion of A, obtained by direct numerical diagonalization, 
with the two analytic limiting expressions and an inter- 
polating ansatz for a chain of L = 100 sites. We will 
now show that the phenomenon of QLDE is intimately 
related to an algebraic scaling behavior of the energy gap 
that is radically different from the one, exponentially de- 
creasing with the size of the chain, exhibited by the fully 
dimerized XX chain. This feature turns out to play a 



Figure 5: (Color online) Behavior of the energy gap Al as 
a function of the length L of the chain for different values of 
A. Diamonds denote the exact numerical data. Continuous 
curves display the approximate expression for the gap (Eq. 
(1331) in the text). 



crucial role at finite temperature, allowing QLDE to be 
more resilient to thermal fluctuations not only at very 
low but even at moderately low temperatures. As al- 
ready pointed out, the single-particle dispersion law has 
the form = Jcosk. The lowest gap is then given by 
J cos (k* ) where k* is the solution of Eq. ([26]) closest to 
7r/2 with k* < 7r/2. For sufficiently large values of L, 
the LHS of Eq. (|26|) has a vertical asymptote arbitrarily 
close to 7r/2. Expanding the function around this point 
we obtain for the desired solution 



k* ~ -- 



2(i-l) (L- 1) 2 A 2 (2- A 2 )" 1 +2(L- 1) 



Then, at leading order in inverse powers of the length of 
the chain L, the lowest energy gap reads 



J 



2(i-l) (L-1) 2 ^ + 2(L-1) 



(33) 



-L~ 



ir(2-\ 2 



L- 2 + O (L- 3 ) 



We have compared the asymptotic analytic behaviors for 
different A and chains of length up to L — 1000 sites with 
the exact numerical solutions. The results of this compar- 
ison are plotted in Fig. ©, showing excellent agreement 
of the exact data with the algebraic scaling ([33]) . 

The existence of a lowest energy gap slowly falling off 
algebraically with the size of the system is to be com- 
pared with the exponentially fast vanishing of the same 
quantity in the case of the dimerized XX chain. Looking 
at the structure of the quasimomenta in the two cases, we 
immediately realize that the reason for this very different 
behavior lies in the absence of complex quasimomenta in 
the spectrum of the XX model with small end bonds. On 
the one hand, this is an undesired feature, because it is 
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0.002 0.004 0.006 0.008 0.01 



T/J 

Figure 6: (Color online) Teleportation fidelity for different 
values of A as a function of the rescaled temperature T/J for 
a XX chain with small end bonds, supporting QLDE. Here 
the length of the chain is fixed at L = 50 sites. The horizontal 
dashed line f — 2/3 is the maximum attainable fidelity using 
a classical teleportation channel. 



exactly the presence of a complex quasimomentum that 
allows for true LDE in the dimerized XX model. On the 
other hand this very same feature allows for true LDE 
and at the same time is responsible for the exponentially 
fast vanishing of the energy gap. 

Moving to finite temperature, the exact numerical eval- 
uation of the fermionic end-to-end correlation function 
x is obtained using Eq. ©. After having determined 
the value of x, we can again use Eqs. (fTU|) and (fTTj) to 
evaluate both the end-to-end concurrence and the tele- 
portation fidelity. In Fig. [6] we report the behavior of 
the teleportation fidelity / for different values of A as a 
function of the normalized temperature T / J for a chain 
of L = 50 spins. We see that the higher the limiting 
value of the fidelity at vanishing temperature, the faster 
/ falls off below the classical threshold / = 2/3. This 
behavior is obviously due to the fact that a larger zero- 
temperature / corresponds to a smaller energy gap, and 
hence the system is more sensitive to the disruptive effect 
of thermal fluctuations. 

We can compare the results reported in Fig. [5] for 
an XX chain with small end bonds with the ones ob- 
tained for an Heisenberg antiferromagnetic chain with 
two weakly interacting end probes [ll|. The teleporta- 
tion efficiency at finite temperature for both models has 
been investigated at the same fixed length of the chain 
L = 50. From this comparison, one sees that the behavior 
of the teleportation fidelity as a function of the tempera- 
ture in the two models is qualitatively similar. This fact 
is of relevance, because it shows that the QLDE due to 
the presence of weakly interacting probes at the and of a 
uniformly interacting chain is an effect not restricted to a 
particular Hamiltonian model, but can be supported by 
systems endowed with different bulk interactions and dif- 
ferent symmetries. The only important ingredient for the 



realization of a QLDE and a teleportation channel that 
are robust against thermal fluctuations is the availability 
of efficient control on the interactions at the end points 
of the channel. Similar conclusions can be drawn about 
LDE: In this work we have demonstrated that existence 
of true LDE at zero temperature is a phenomenon not 
restricted to a particular model and, in fact, does not 
crucially depend on the mechanism of dimerization. Pre- 
liminary studies indicate that gapped, anisotropic mod- 
els in non critical regimes can sustain LDE even in the 
absence of a clear pattern of dimerization. 

Finally, it is worth noticing that even if the Heisen- 
berg model with small end bonds analyzed in Ref. [llj 
sustains a teleportation fidelity that remains above the 
classical threshold for larger values of the end coupling 
parameter A, the XX channel with small end bonds as- 
sures teleportation with nonclassical fidelities for larger 
values of the temperature. Depending on the realistic ex- 
perimental settings, if the temperature can be controlled 
down to sufficiently low values, it is more convenient to 
engineer a channel based on the Heisenberg-like Hamil- 
tonian proposed in Ref. Otherwise, if control on the 
external temperature cannot be implemented with high 
precision down to sufficiently low temperatures, it may 
then turn out more convenient to engineer XX chains of 
the type discussed in the present paper. 



V. CONCLUSIONS AND OUTLOOK 

In conclusion, we have studied models of open quan- 
tum spin chains endowed with XX-like Hamiltonians 
with nearest-neighbor interactions. We have discussed 
two types of models, according to different engineerings 
of the spin-spin couplings. In the first case considered, a 
chain with bonds of alternating strengths, we have shown 
how the exact solution for the ground state implies the 
existence of long-distance entanglement between the end 
spins of the chain, independent of the size of the sys- 
tem and asymptotically close to unity (maximal entan- 
glement) in the limit of exact dimerization. This system 
is therefore perfectly suited for bona fide long-distance 
quantum teleportation with ideal fidelity at zero tem- 
perature. However, the limiting maximal values of the 
fidelity are obtained for an exponentially small energy 
gap above the ground state. Therefore, this system is de 
facto useless for efficient quantum teleportation at finite 
temperature. We have then moved on to discuss another 
class of XX open spin chains with uniform bulk interac- 
tions and small end bonds. In this case, we have shown 
that for sufficiently small values of the end couplings, the 
ground state of the system supports a quasi long-distance 
entanglement between the end spins of the chain, asymp- 
totically close to unity (maximal entanglement) in the 
limit of vanishing coupling, but slowly decreasing as the 
length of the chain is increased. An interesting feature of 
this model is that the lowest energy gap above the ground 
state vanishes only algebraically, as the first power of the 
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inverse of the size of the system. Therefore, in principle, 
it can be exploited as a quantum channel for teleporta- 
tion with nonclassical fidelity at finite temperature, both 
very low and moderately low. 

Considering future research along these lines of investi- 
gation, it will be interesting to consider practical schemes 
for the realization of this kind of spin chains in highly 
controllable situations, for instance resorting to ultracold 
atomic mixtures in optical lattices. Another interesting 
open problem worth further study is to assess the exis- 
tence and the possible location of a crossover between 
true long-distance and prima facie quasi long-distance 



entanglement behaviors. 
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